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METHOD OF COMPUTING THE PITCH NAMES 
OF NOTES IN MIDI-LIKE MUSIC 
REPRESENTATIONS 

Field of the invention 

This invention relates to the problem of constructing a reliable pitch spelling algorithm — 
that is, an algorithm that reliably computes the correct pitch names (e.g., CJJ4, Bb5 etc.) 
of the notes in a passage of tonal music, when given only the onset-time, MIDI note 
number (The MIDI Manufacturers' Association, 1996, p. 10) and possibly the duration of 
each note in the passage. 

There are good practical and scientific reasons for attempting to develop a reliable pitch 
spelling algorithm. First, until such an algorithm is devised, it will be impossible to con- 
struct a reliable MIDI-to-notation transcription algorithm — that is, an algorithm that reli- 
ably computes a correctly notated score of a passage when given only a MIDI file of the pas- 
sage as input. Commercial music notation programs (e.g., Sibelius (www.sibelius.com), 
Coda Finale (www. codamusic . com) and Nightingale (www.ngale.com)) typically use 
MIDI-to-notation transcription algorithms to allow the user to generate a notated score 
from a MIDI file encoding a performance of the passage to be notated. However, the 
MIDI-to-notation transcription algorithms that are currently used in commercial music 
notation programs are crude and unreliable. Also, existing audio transcription systems 
generate not notated scores but MIDI-like representations as output (see, for example, 
Davy and Godsill, 2003; Plumbley et al : 2002; Walmsley, 2000). So if one wishes to 
produce. a notated score from a digital audio recording, one typically needs a MIDI-to- 
notation transcription algorithm (incorporating a pitch spelling algorithm) in addition to 
an audio transcription system. 

Knowing the letter-names of the pitch events in a passage is also indisp.ensible in music 
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information retrieval and musical pattern discovery (Meredith et al, 2002). For example, 
the occurrence of a motive on a different degree of a scale (e.g., C-D-E-C restated as 
E-F-G-E) might be perceptually significant even if the corresponding chromatic intervals 
in the patterns differ. Such matches can be found using fast, exact-matching algorithms if 
the pitch names of the notes are encoded, but exact-matching algorithms cannot be used 
to find such matches if the pitches are represented using just MIDI note numbers. If a 
reliable pitch spelling algorithm existed, it could be used to compute the pitch names of 
the notes in the tens of thousands of MIDI files of works that are freely available online, 
allowing these files to be searched more effectively by a music information retrieval (MIR) 
system. 

In the vast majority of cases, the correct pitch name for a note in a passage of tonal 
music can be determined by considering the roles that the note is perceived to play in the 
perceived harmonic structure and voice-leading structure of the passage. For example, 
when played in isolation in an equal-tempered tuning system, the first soprano note in 
Figure la would sound the same as the first soprano note in Figure lb. However, in 
Figure la, this note is spelt as a Gj}4 because it is perceived to function as a leading note 
in A minor; whereas in Figure lb, the first soprano note is spelt as an Ab4 because it 
functions as a submediant in C minor. Similarly, the first alto note in Figure lb would 
sound the same as the first alto note in Figure lc in an equal-tempered tuning system. 
However, in Figure lb the first alto note is spelt as an Fij4 because it functions in this 
context as a subdominant in C minor; whereas, in Figure lc, the first alto note functions 
as a leading note in Ftf minor so it is spelt as an EJJ4. 

Nevertheless, it is not always easy to determine the correct pitch name of a note 
by considering the harmonic structure and voice-leading structure of its context. For 
example, as Piston (1978, p. 390) observes, the tenor Eb4 in the third and fourth bars of 
Figure 2 should be spelt as a D}14 if one perceives the harmonic progression here to be 
+II 2 — I as shown. But spelling the soprano Eb5 in the fourth bar as Dft5 would result in 
a strange melodic line. 
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Such cases where it is difficult to determine the correct pitch name of a note in a tonal 
work are relatively rare— particularly in Western tonal music of the so-called 'common 
practice 5 period (roughly the 18th and 19th centuries). In the vast majority of cases, those 
who study and perform Western tonal music agree about how a note should be spelt in 
: a given tonal context. Therefore a pitch spelling algorithm can be evaluated objectively 
by running it on tonal works and comparing the pitch names it predicts with those of the 
corresponding notes in authoritative published editions of scores of the works. However, 
this can only be done accurately and quickly if one has access to encodings of these 

, -t . , . , . . •„ j-T r ^ filoc +V.o+ ^o-n Ko ^mnflrpH p 1 1 torn pi ,t i r: a,] l.V 

aUUQOIiCaxrve t^oies m b-LLC iuiiii wi wniyuuoj. ~~ ~ r- - " 

with the pitch spelling algorithm's output. 

Related Art 

In this section, the performance of three prior pitch spelling algorithms is compared on 
a single test corpus of works containing 41544 notes and consisting of all 48 pieces in 
the first book of J. S. Bach's Das Wohltemperirte Klavier. The algorithms compared are 
those of Cambouropoulos (1996, 1998, 2000, 2001, 2002), Longuet-Higgins (1976, 1987, 
1993) and Temperley (1997, 2001). 

When these algorithms were run on the test corpus, Cambouropoulos's algorithm 
made 2599 mistakes, Longuet-Higgins's algorithm made 265 mistakes and Temperley's 
algorithm made 122 mistakes. 

The test corpus The test corpus used in the comparison consists of representations 
of the scores of all 24 Preludes and all 24 Fugues in the first book of J. S. Bach's Das 
Wohltemperirte Klavier encoded in the author's OPND format. 1 Each OPND representa- 
tion is a set of triples, (£, n, d), each triple giving the onset time, t, the pitch name, n, and 
^OPND stands for "onset, pitch-name, duration". 
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the duration, of a single note (or sequence of tied notes) in the score. 2 The onset time 
and duration of each note are expressed as integer multiples of the largest common divisor 
of all the notated onset times and note durations in the score. For example, Figure 3b 
gives the OPND representation of the score in Figure 3a. Note that within each pitch 
name in an element in an OPND representation, each flat symbol is represented by an 
'F 5 character and each sharp- symbol is represented by an 'S' character (double-sharps are 
denoted by two C S' characters, so, for example, F>:4 is denoted by "FSS4".) 

The test corpus was derived by automatic conversion from Hewlett's (1997) MuseData 
encodings of Bach ; s Das Wohiiemperirte Kiavierr The MuseData encodings contained 
minor errors which were corrected in the OPND representations. 

Also, Temperley's algorithm cannot deal with situations in which two or more notes 
with the same pitch begin at the same time. Thus, wherever two or more notes, with the 
same pitch p began simultaneously at time £, all notes with pitch p and onset time t were 
removed from the OPND file except the one with the longest duration! This resulted in 
a test corpus containing 41544 notes. 4 £ 

The "piano-roll" or MIDI-like input representations accepted by the algorithms com- 
pared in this study were derived automatically from the OPND representations of the 
pieces in the test corpus. 

Longuet-Higgins's algorithm Pitch spelling is one of the tasks performed by Longuet- 
Higgins's (1976, 1987, 1993) music. p program. Longuet-Higgins has published the full 
POP-11 source code of this program (Longuet-Higgins, 1987, pp. 120-126; Longuet- 
Higgins, 1993, pp. 486-492). Just the pitch spelling portion of music.p was translated 
directly into Lisp in order to make it easier to perform the comparison. This was possi- 
ble because the pitch spelling portion of music.p operates independently of the rhythmic 

2 The voice of each note is also given in the files used in the comparison, but this voice information is 
not used by any of the algorithms. 

3 Available online at http://ww.musedata.org/encodings/bacii/bg/keybd/. 

4 The complete test corpus is available online at http: //www. tit anmusic - com/ data. html. 
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part. 5 

The input to music. p must be in the form of a list of triples, (p, ton, *off), eac k triple 
giving the "keyboard position" p together with the onset time t on and the offset time £ 0 ff 
in centiseconds of each note. The keyboard position p is simply an integer indicating the 
key that would have to be pressed on a normal piano keyboard in order to perform the 
note, with C\\3 mapping onto 0, Cp and Db3 mapping onto 1, CH4 mapping onto 12 and 
so on (i.e., p = MIDI NOTE NUMBER - 48). 

The algorithm then computes a value of "sharpness" q for each note in the input 
(Longuet-Higgiiis, 1987, p. 111). The sharpness of a pitch name indicates the position of 
the pitch name on the line of fifths (Temperley, 2001, p. 117) and is therefore essentially the 
same as Temperley's (2001, p. 118) concept of "tonal pitch class". In Longuet-Higgins 5 s 
algorithm, if <?i and q 2 are the sharpnesses of two notes then the interval between the 
notes is defined to have "degree" 6q = q 2 - ft. If \Sq\ < 6, the interval is defined to be 
"diatonic"; if \6q\ > 6, it is defined to be "chromatic"; and if \Sq\ = 6, it is defined to be 
"diabolic" (Longuet-Higgins, 1987, p. 112). Longuet-Higgins's algorithm attempts to spell 
notes so that the degree between each note and the tonic at the point at which the note 
occurs is not chromatic (Longuet-Higgins, 1987, p. 113). The algorithm also incorporates 
various rules for dealing with chromatic passages (Longuet-Higgins, 1987, pp. 113-114). 

Pitch spelling algorit hms sometimes spell complete pieces in a different key from that 
in which they are notated in the original score. For example, Longuet-Higgins's algorithm 
spells the Prelude in Gfi minor from Book 1 of Das Wohltemperirte Klavier (BWV 863) 
in Ab minor. This should not be considered an error. What matters is whether or not the 
pitch interval names (e.g., 'rising major third 5 , 'falling augmented fourth 5 , etc.) between 
corresponding pairs of notes are the same in the computed spelling as they are in the 
published score. The author's implementation of Longuet-Higgins's algorithm therefore 
actually generates three alternative spellings for each piece: 

5 The Lisp implementation of the pitch spelling portion of music.p is available online at 
http : //www . titanmusic. com/ soft ware . html. 
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1. the spelling S for the whole input passage computed by the algorithm; 

2. the spelling that results when S is transposed by a rising diminished second; and 

3. the spelling that results when S is transposed by a falling diminished second. 

Each of these three computed spellings is then compared with the pitch names in the 
original score and the number of errors made by the program is taken to be the number 
of errors in the best of the three alternative spellings. 

When the author's implementation of Longuet-Higgins's algorithm was run on the test 
corpus described above, it made only 265 errors-that is, it predicted the correct pitch 
name for 99.36% of the notes. 

Cambouropoulos's algorithm Unlike Longuet-Higgins, Cambouropoulos has not 
published an implementation of his pitch spelling algorithm, nor was. he able to sup- 
ply his own implementation when it was requested. A new implementation of his method 
was therefore made, based on his published descriptions of it (Cambouropoulos, 1996, 

1998, 2000, 2001, 2002). 6 

■ , . - -. . -. .j. , «. +Uo 5mvn+ T-onr<^<3ftr)t.ation into a 

Cambouropoulos's metnod. mvoives nrst cuxivcumg ^ — <Lj -- 

sequence of chromas or pitch classes in which the chromas are in the order in which 
they occur in the music (the chromas of notes that occur simultaneously being ordered 
arbitrarily). The chroma of a note can be computed from its MIDI note number using 
the formula 

CHROMA = MIDI NOTE NUMBER mod 12. 

The term chroma has been used in this sense since around 1950 (Bachem, 1950; Burns 
and Ward, 1982, pp. 246, 262-264; Cross et al, 1991, pp. 212, 223-224; Deutsch, 1982, 
p. 272; Deutsch, 1999, p. 350; Shepard, 1964; Shepard, 1965; Shepard, 1982, p. 352; 
Dowling, 1991, p. 35; Ward and Burns, 1982, p. 432-433 ). The term chroma is essentially 

6 T he Lisp code for the author's implementation of Cambouropoulos's pitch spelling algorithm is avail- 
able online at http://ww.titanmusic.com/softwar-e.html. 
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"synonymous with the term pitch class as used in atonal theory (Babbitt, 1960; Babbitt, 
1965; Forte, 1973; Morris, 1987; Rahn, 1980). 

Having derived an ordered set of chromas from the input, Cambouropoulos's algorithm 
then processes the music a window at a time, each window containing a fixed number of 
notes (set to a value between 9 and 15). Each window is positioned so that the first third 
of the window overlaps the last third of the previous window. Cambouropoulos allows 
'white note' chromas (i.e., 0, 2, 4, 5, 7, 9 and 11) to be spelt in three different ways 
(e.g., chroma 0 can be spelt as BJt, Q or Dbb) and 'black note' chromas to be spelt in 
two different ways (e.g., chroma 6 can be speit as Fjj or Gb). Given these restricted sets 
of possible pitch names for each chroma, the algorithm computes all possible spellings 
for each window. A penalty score is then computed for each of these, possible window 
spellings. The penalty score for a given window spelling is found by computing a penalty 
value for the interval between each pair of notes in the window and summing these penalty 
values. A given interval in a particular window spelling is penalised more heavily if it is 
an interval that occurs less frequently in the major and minor scales. An interval is also 
penalised if either of the pitch names forming the interval is a double-sharp or a double- 
flat. For each window, the algorithm chooses the spelling that has the lowest penalty 



score. 



When the author's implementation of Cambouropoulos's method was run on the test 
corpus, it made 2599 mistakes— that is, it predicted the correct pitch name for only 93.74% 
of the notes. When Cambouropoulos ran his own implementation of his method on the 
test corpus, he found that it made even more errors than the author's implementation. 
This may be due to the fact that the new implementation -generates three alternative 
transpositions of the computed spelling and chooses the one that results in the least 
number of errors. 



Temperley's algorithm Temperley's (1997, 2001) pitch spelling algorithm is imple- 
mented in his harmony program which forms one component of his and Sleator's Melisma 
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system. 7 The input to the harmony program must be in the form of a "note-list" (Tem- 
perley, 2001, pp. 9-12) giving the MIDI note number of each note together with its onset 
time and duration in rnilliseconds. The harmony program also requires a specification 
of the metrical structure of the input passage which can be generated by first running 
the note-list through Temperley and Sleator's meter program (another component of the 
Melisma system). 

Temperley's (2001, pp. 115-136) pitch spelling algorithm searches for the spelling 
that best satisfies three "preference rules". The first of these rules stipulates that the 
algorithm should "prefer to label nearby events so that they axe close together on the 
fine of fifths" (Temperley, 2001, p. 125). This rule bears some resemblance to the basic 
principle underlying Longuet-Higgins's algorithm (see above). The second rule expresses 
the principle that if two tones are separated by a semitone and the first tone is distant from 
the key centre, then the interval between them should preferably be spelt as a diatonic 
semitone rather than a chromatic one (Temperley, 2001, p. 129). The third preference rule 
'underlying Temperley's algorithm steers the algorithm towards spelling the notes so that 
a "good harmonic representation" results (Temperley, 2001, p. 131), a "good harmonic 
representation" being one that allows Temperley's harmony program to generate a correct 
harmonic analysis of the passage. 
- Because the output of the harmony program depends on tempo, each OPND repre- 
sentation in the test corpus had to be supplemented with a file giving a tempo map for 
the piece represented. These two files were then automatically converted into a note-fist 
in the form required by Temperley's programs. The output of the harmony program 
was then automatically converted into a format which would allow automatic comparison 
with the test corpus files. Again, for each piece in the test corpus, three transpositions of 
. the spelling generated by Temperley's program were compared with the original spelling 
and the one with the least number of errors was selected. When Temperley's harmony 
'The complete Melisma system together with documentation is available online at 
http : / /www . link. cs . emu . edu/nmsic-analysis/ . 
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program was run on the test corpus, it made only 122 mistakes — that is, it predicted the 
correct pitch name for 99.71% of the notes. 

Summary of the invention 

The invention described here consists of an algorithmic method called pslS that reliably 
computes the correct pitch names (e.g., C$4, Bb5 etc.) of the notes in a passage of tonal 
music, when given only the onset-time and MIDI note number of each note in the passage. 

The pslS algorithm has been shown to be more reliable than previous algorithms, 
correctly predicting the pitch names of 99.81% of the notes in a test corpus containing 
41544 notes and consisting of all the pieces in the first book of J. S. Bach's Das Wohltem- 
perirte Klavier (Le:, pslS incorrectly predicted the pitch names of only 81 notes in this 
test corpus). 

Three previous algorithms (those of Cambouropoulos (1996, 1998, 2000, 2001, 2002), 
Longuet-Higgins (1976, 1987, 1993) and Temperley (1997, 2001)) were run on the same 
corpus of 41544 notes. On this corpus, Cambouropoulos's algorithm made 2599 mistakes, 
Longuet-Higgins's algorithm made 265 mistakes and Temperley's algorithm made 122 

mistakes. As pslS made only 81 mistakes on the same corpus, this provides evidence in 

< 

support of the claim that pslS is more reliable than previous algorithms that attempt to 
perform the same task. 

The pslS algorithm is best understood to be in two parts, Part I and Part II. Part I 
consists of the following steps: 

1. computing for each pitch class 0 < c < 11 and each note n in the input, the pitch 
letter name S{c y n) E {A, B y C, D, E, F, G} that n would have if c were the tonic at 
the point in the piece where n occurs (assuming that the notes are spelt as they are 
in the harmonic chromatic scale on c); 

2. computing for each note n in the input and each pitch class 0 < c < 11, a value 
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CNT(c, n) giving the number of times that c occurs within a context surrounding n 
that includes n, some specified number K pie of notes immediately preceding n and 
some specified number X post of notes immediately following n; 

3. computing for each note n and each letter name I, the set of chromas C(n,l) = 
{c | S(c,n) = 1} (that is, the set of tonic chromas that would lead to n having the 
letter name Z); 

4. computing N(l, n) = £ CJVT(c, n) for each note n and each pitch letter name I; 

cEC(n,l) 

5. computing for each note n, the letter name I for which N(l, n) is a maximum. 

Part II of the algorithm corrects those instances in the output of Part I where a 
neighbour note or passing note is erroneously predicted to have the same letter name as 
either the note preceding it or the note following it. Part II of pslS 

1. lowers the letter name of every lower neighbour note for which the. letter name 
predicted by Part I is the same as that of the preceding note; 

2. raises the letter name of every upper neighbour note for which the letter name 
predicted by Part I is the same as that of the preceding note; 

3. lowers the letter name of every descending passing note for which the letter name 
predicted by Part I is the same as that of the preceding note; 

4. raises the letter name of every descending passing note for which the letter name 
predicted by Part I is the same as that of the following note; 

5. lowers the letter name of every ascending passing note for which the letter name 
predicted by Part I is the same as that of the following note; 

6. raises the letter name of every ascending passing note for which the letter name 
predicted by Part I is the same as that of the preceding note. 
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Detailed description of preferred implementations 

The invention consists of an algorithm, called psl3, that takes as input a representation of 
a musical passage (or work or set of works) in the form of a set I of triples (t,p c , d) , each of 
these triples giving the onset time t, the chromatic pitch p c and the duration d of a single 
note or sequence of tied notes. The chromatic pitch of a note is an integer indicating the 
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interval in semitones from A*0 to the pitch of the note (i.e., p c = MIDI NOTE NUMBER— 
21). The set I may be trivially derived from a MIDI file representation of the musical 
passage. 

If c, = ft, Pcii *) and e, = ft, «V> and «, e, € I then e, is denned to be less than 
ej, denoted by e, < e j} if and only if U < h or (U = tj Ap c>i < Pe j) or (U = t s Ap C) i = 
p cJ A(U< dj). Also, e« < e,- if and only if e. : < e, or e, = e,. 

The first step in ps^5 is to sort the set J to give an ordered set 

dui» (D 

containing all and only the elements of J, sorted into increasing order so that j > i 

<*i,Pci, di) < (*i,ftj> d i> for aU di > ' ( *'"' PcJ '' di) 6 J - 

The second step in psi3 is to compute the ordered set 



C = <ci,c 2 ,...c fc , ...C|j|) 



(2) 



where Cfc = pcfe mod 12, PcM being the chromatic pitch of the fcth element of J. c k is the 

chroma of p c ,k- 

If A is an ordered set of elements, 

A= (a 1 , a 2 , ■ ■ ■ a k , ■ ■ ■ a\A\) 

then let A\j] denote the (j + l)th element of A (e.g., A[0] = a 1; A[l] = a 2 ). Also, let 
A\j, k] denote the ordered set that contains all the elements of A from A\j) to A[k - 1], 
inclusive (e.g., A[l, 4] = {a 2 ,a 3 ,a 4 )). 

In addition to the set I, pslS takes as input two numerical parameters, caUed the 
precoma denoted by K pie> and the postcard**, denoted by iW The precontext must 
be an integer greater than or equal to 0 and the postcontext mnst be an integer greater 
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than 0. The third step in pslS is to compute the ordered set 

$ = (4>(C[0)),4>(C[1}), . . . 4>{C[k]), . . . <KC[\C\ - 1])) (3) 

where 

ft C [k]) = {CNT(0, k), CNTQ-, k), . . . CNT{11, k)) (4) 

and CNT(c,k) returns the number of times that the chroma c occurs in the ordered set 
C[max({0 ) fc-K pre }),min({|C!,fc + K P ost})] (the functions max(B) and min(S) return 
the greatest and least values, respectively, in the set B). 

Every pitch name has three parts: a letter name which must be a member of 
the set {A,B,C,D,E,F,G}; an inflection which must be a member of the infinite set 
{. xx, ftx, x, Jj, ||, b, bb, bbb, bbbb, . . .}; and an octave number which must be an integer. By 
convention, if the inflection of a pitch name is equal to t) it may be omitted, for example, 
Q4 may be written C4. The octave number of 'middle C is 4 and the octave number of 
any other C is one greater than the next C below it on the staff and one less than the 
next C above it on the staff. The octave number of any pitch name whose letter name is 
not C is the same as that of the nearest C below it on the staff. Thus Bx3 sounds one 
semitone higher than C4 and Cb4 has the same sounding pitch as Bt|3. 

If iV is a pitch name with letter name l(N) and octave number o(N), then the morph 
of N, denoted by m(N), is given by the following table 



l(N) 


A B C D E F G 


m(N) 


0 1 2 3 4 5 6 



The morphetic octave of TV, denoted by o m (N), is given by 



o(N), if m = 0 or rrc = 1, 
o(N) — 1, otherwise. 
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The morphetic pitch of N, denoted by p m O")> is given by 

p m (N) = m(N) + 7o m (JV). (6) 

The fourth step in V sl3 is to compute the spelling table, S. This is done using the 
algorithm expressed in pseudo-code in Figure 4. In this pseudo-code, block structure is 
indicated by indentation and the symbol «-» is used to denote assignment of a value to 
a variable (i.e., the expression «x <- y" means "the value of variable x becomes equal to 
the value of if). The symbol "©" is used to denote concatenation of ordered sets. Thus, 
if A and B are two ordered sets such that A = ~(a u a 2 , . . . a ]A \) and B = (h, b 2 , . . . b m ), 
then 

A®B = (a 1 ,a 2 ,..-aiA\,bi,b2,.--b\B\)- 

Also, A®() = ()®A = A. 

Having computed the spelling table S, the algorithm goes on to compute the relative 
morph list, R. This is computed using the algorithm in Figure 5. If A is an ordered set 
of ordered sets, then A[i\\j] denotes the (j + l)th element of the (i + l)th eiement of A. 
The expression S[k)\i] in Une 6 of Figure 5 therefore denotes the (i + l)th element of the 
(k + l)th element of the spelling table S. If L is an ordered set and i is the value of at 
least one element of L, then the function POS(i,L) called in line 10 of Figure 5 returns 
the least value of k for which L[k] = i. 

The next step in pslS is to compute the initial morph m init which is given by 

Qu*[C[0]] ( ? ) 

where Q Snit = (0,1,1,2,2,3,4,4,5,5,6, 6). 

Next, psl3 computes the morph list, M, which satisfies the following condition 



(\M\ = \C\) A (M[i] = (i?[«] + minit) mod 7 for all 0 < i < \C\). 



(8) 
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Next, the ordered set O is computed which satisfies the following condition 

(|0| = |C|) A (0[i] = <J[i][0], C[i], M[i]> for all 0 < i < \C\). (9) 

Next, an ordered set # called the chord list is computed using the algorithm in Fig- 
ure 6. 

Ia the next three steps, the algorithm ensures that neighbour note and passing note 
figures are spelt correctly. Figure 7 illustrates the six types of error corrected by the 
algorithm. Figure 7a shows a lower neighbour note figure that is incorrectly spelt because 
the neighbour note (the middle note) has the same morphetic pitch as the two flanking 
notes. In this case, the morphetic pitch of the neighbour note is one greater than it should 
be. Figure 7b shows a similar error in the case of an upper neighbour note figure. In this 
case, the morphetic pitch of the neighbour note is one less than it should be. The passing 
note figures in Figure 7c and d are incorrect because the morphetic pitch of the passing 
note (the middle note) is one greater than it should be. Finally, the passing note figures 
in Figure 7e and f are incorrect because the morphetic pitch of the passing note is one 

less than it should be. 

pslS first corrects errors like the ones in Figure 7a and b using the algorithm in 
Figure 8. If A is an ordered set of ordered sets then the expression A[i][j, fc] denotes the 
ordered set (A\i)\j], A[i\[j + 1], • ■ ■ A[i\[k *- l]>. Thus, the expression H[i + 2)[k)[l,3) in 
fine 4 of Figure 8 denotes the ordered set (H[i + 2}{k}[l] ) H[i + 2}[k}[2)). In Figure 7a and 
b, the neighbour note is one semitone away from the flanking notes. The algorithm in 
Figure 8 also corrects instances where the neighbour note is 2 semitones above or below 
the flanking notes but has the same morphetic pitch as the flanking notes. 

Next, pslS corrects errors like the ones in Figure 7c and e using the algorithm in 
Figure 9. Then it corrects errors like the ones in Figure 7d and f using the algorithm in 
Figure 10. In Figure 7c, d, e and f, the interval between the flanking notes is a minor third. 
The algorithms in Figures 9 and 10 also correct instances where the interval between the 
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flanking notes is a major third. 

Having corrected neighbour note and passing note errors, psl 3 then computes a new 
morph list M' using the algorithm in Figure 11, 

Having computed M', it is now possible to compute a morphetic pitch for each note. 
This can' be done using the algorithm in Figure 12 which computes the ordered set of 
morphetic pitches, P. 

Finally, from J and P, psl 3 computes an OPND representation Z which satisfies the 

following condition 

{\Z\ = \J\) A (Z\i] = (J[t][0], PPN({J[i\[l), P[i\)), JT«J[2]> for aU 0 < i < \Z\). (10) 

The function PPNdpcPm}) returns the unique pitch name whose chromatic pitch is p c 
and whose morphetic pitch is p m . This function can be computed using the algorithm in 
Figure 13. In this algorithm, anything in double inverted commas ("?" ) is a string— that 
is, an ordered set of characters. If A and B are strings such that A = "abode f" and 
B = "ghijkl", then the concatenation of A onto B, denoted by A®B, is " abode f g hi jkF , 
In the pitch names generated by PPN, the sharp signs are represented by <s' characters and 
the flat signs are represented by 'f ' characters. A double sharp is represented by the string 
"ss". For example, the pitch name Cx4 is represented in PPN by the string "Css4". The 
empty string is denoted by and the function STR(n) called in line 23 returns a string 
representation of the number n. For example, STR(-5) = " - 5" , STR(105.6) = "105.6" . 

Results of running ps 1 3 on the test corpus As explained above, psl 3 takes as input 
a set / of triples, (t,p c ,d), each one giving the onset time, chromatic pitch and duration 
of each note. In addition, psl 3 requires two numerical parameters, the precontext, K pTey 
and the postcontext, K post - ' . ■ * 

In order to determine the values of K pre and if post that give the best results, psl 3 
was run on the test corpus 2500 times, each time using a different pair of values 
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(K pte , tf pos t> chosen from the set {(K pie , K post ) | 1 < K pt *,K post < 50}. For each pair of 
values (K pie , K post ), the number of errors made by pslS on the test corpus was recorded. 

pslS made fewer than 122 mistakes (i.e., performed better than Temperley's algorithm) 
on the test corpus for 2004 of the 2500 (K pi&) K^) pairs tested (i.e., 80.160% of the 
(K pte , impost) pairs). 

pslS performed best on the test corpus when K pTS was set to 33 and i<r post was set 
to either 23 or 25. With these parameter values, pslS made only 81 errors on the test 
corpus— that is, it correctly predicted the pitch names of 99.805% of the notes in the test 
corpus. 

The mean number of errors made by pslS over all 2500 (K pie , K post ) pairs was 109.082 
(i.e., 99.737% of the notes were correctly spelt on average over all 2500 (K pie , K post ) pairs). 
This average value is better than the result obtained by Temperley's algorithm for this 
test corpus. 

The worst result was obtained when both K pie and K" pos t were set to 1. In this case, 
pslS made 1117 errors (97.311% correct). However, provided K pie is greater than about 
14 and Kpost is greater than about 21, pslS predicts the correct pitch name for over 
99.75% of the notes in the test corpus. 

A Lisp implementation of psl3 The Lisp implementation of pslS given below as- 
sumes that the input representation I is represented as a list of sublists, each sublist 
taking the form (t p c d) where t, p c and d are the onset time, chromatic pitch and du- 
ration, respectively, of a single note. For example, the passage in Figure 3a would be 
represented by the list 

((2 41 2) (4 43 2) (6 44 2) (8 46 2) (10 43 2) (12 44 1) (13 41 1) (14 40 1) 

(15 41 1) (16 49 4) (20 46 4) (24 48 5) (26 36 2) (28 38 2) (29 46 1) 

(30 39 2) (30 44 1) (31 43 1) (32 41 2). (32 46 1) (33 44 1) (34 38 2) 

(34 43 1) (35 41 1) (36 39 1) (36 43 2) (37 36 1) (38 35 1) (38 51 3) 

(39.36 1) (40 44 4) (41 50 1) (42 48 1) (43 50 1) (44 41 4) (44 50 1) 

(45 48 1) (46 47 1) (47 48 1)) 

Similarly, the output of the Lisp implementation of pslS given below is represented as 
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a list of sublists. For example, the output of this implementation of pslS for the passage 
given in Figure 3a is 

((2 "Dn4" 2) (4 "En4" 2) (6 "Fn4" 2) (8 "Gn4" 2) (10 "En4" 2) (12 "Fn4" 1) 

(13 "Dn4" 1) (14 "Cs4" 1) (15 "Dn4" 1) (16 "Bf4" 4) (20 "Gn4" 4) (24 "An4" 5) 

(26 "An3" 2) (28 "Bn3" 2) (29 "Gn4" 1) (30 "Cn4" 2) (30 "Fn4" 1) (31 "En4" 1) 

,r>~ »n.«» /-oo »G"4» 1 "> '3^ "Vr,4." 1 ~> f34 "Br>3" 2") (34 "En4" 1) (35 "Dn4" 1) 

(36 "Cn4" 1) (36 "En4" 2) (37 "An3" 1) (38 "Gs3" 1) (38 "Cn5" 3) (39 "An3" 1) 

(40 "Fn4" 4) (41 "Bn4" 1) (42 "An4" 1) (43 "Bn4" 1) (44 M Dn4 M 4) (44 "Bn4" 1) 
(45 "An4" 1) (46 "Gs4" 1) (47 "An4" 1)) 



Mere, then, is the Lisp code for an impIej-Qexxtation of psl 3: 

(defun psl3 (feoptional (input-filename (choose-f ile-dialog) ) 

(pre-context 33) 
(post-context 23)) 
(let* ((sorted-input-representation 
(remove -duplicates 
(sort (with-open-f ile (input-file 

input - f i 1 ename ) 
(read input-file)) 
# ' vector-less-than) 
:test #'equalp)) 
(onset-list (mapcar #' first sorted-input-representation)) 
( chr omat i c-p it ch-1 i st 

(mapcar #' second sorted-input-representation)) 
(chroma-list 

(mapcar #' chromatic -pitch-chroma chromatic-pitch-list)) 
(n (list-length chroma-list)) 
( chroma- ve ct or-1 is t 
(do* ((cvl nil) 

(i 0 (1+ i))) 
((= i n) 
cvl) 
(setf cvl 

(append cvl 

(list 
(do* ((context 

(subseq chroma-list 

(max 0 (- i pre-context)) 
(min n (+ i post-context)))) 
(cv (list 0 0 0 0 0 0 0 0 0 0 0 0)) 
(c 0 (+ 1 c))) 
((= c 12) 
cv) 

(setf (elt cv c) 
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(count c context)))))))) 
(chromamorph-table (list 01122334556 6)) 
(spelling-table 
(do* ( (first -morph nil nil) 
(spelling nil nil) 
(spelling2 nil nil) 
(st nil) 
(c 0 (1+ c))) 
((= c 12) 
st) 

(setf spelling 

(mapcar #' (lambda (chroma-in-chroma-list) 
(elt chromamorph-table 

(mod (- chroma- in- chroma- list c) 12))) 
chroma- list) ) 
(setf first-morph (first spelling)) 
(setf spelling2 

(mapcar # J (lambda (morph- in- spelling) 

(mod (- morph- in-spelling first-morph) 7)) 
spelling) ) 

(setf st (append st (list spelling2) ) ) ) ) 
(relative-morph-list 
(do ((morph- vector (list 0 0 0 0 0 0 0) 
(list 0 0 0 0 0 0 0)) 

(rml nil) 
(i 0 (1+ i)) 

(morphs-f or-this- chroma nil 

nil) 

) 

((= i n) 
rml) 

(setf morphs-f or-this-chroma 

(mapcar #' (lambda (spelling) 
(elt spelling i)) 
spelling-table) ) 

(setf rml 

(do ((prev-score nil nil) 
(j 0 (1+ j))) 
((= j 12) 
, ; (pprint morph-vector) 
(append rml 

(list (position 

(apply #'max morph-vector) 
morph-vector) ) ) ) 

(setf prev-score 

(elt morph-vector 

(elt morphs-f or-this-chroma j))) 
(setf (elt morph-vector 
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(elt morphs-f or-this- chroma j)) 
(+ prev-score 

(elt (elt chroma-vector-list i) j))))))) 
(initial-morph (elt '(0 1122344556 6) 

(mod (first chromatic -pit ch-list). 12))) 
(morph-list (mapcar #' (lambda (relative-morph) 

(mod (+ relative-morph initial-morph) 7)) 
relative-morph-list) ) 
(ocm (mapcar # 3 list onset -list chroma-list morph-list)) 
(ocm-chord-list (do* ((cl (list (list (first ocm)))) 

(i 1 (1+ i))) 
((= i n) 
cl) 

(if (= (first (elt ocm i)) 

* /■* _ 

(setf (first (last cl)) 

(append (first (last cl)) 

(list (elt ocm i)))) 

(setf cl 

(append cl 

(list (list (elt ocm i) ))))))) 

(number-of -chords (list-length ocm-chord-list)) 

; neighbour notes 

(ocm-chord-list 

(do* ((i 0 (1+ i))) 

(0= i (- number-of -chords 2)) 
ocm-chord-list) 
(dolist (notel (elt ocm-chord-list i)) 
(if (member (car notel) 

(mapcar #'cdr (elt ocm-chord-list (+ i 2))) 

ttest #>equalp) 
(dolist (note2 (elt ocm-chord-list (1+ i))) 
(if (= (third note2) 
(third notel)) 
(progn 

(if (member (mod (- (second note2) (second notel)) 
12) 
'(1 2)) 
(setf (third note2) 

(mod (+ 1 (third note2)) 7))) 
(if (member (mod (- (second notel) (second note2)) 
12) 
'(1 2)) 
(setf (third note2) 

(mod (- (third note2) 1) 7)))))))))) 

; downward passing notes 
(ocm-chord-list " 
(do* ((i 0 (1+ i))) 
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((- i (- number-of -chords 2)) 
o cm- chord-list ) 
(dolist (notel (elt ocm-chord-list i)) 

(dolist (note3 (elt ocm-chord-list (+ i 2))) 

(if (- (third note3) (mod (- (third notel) 2) 7)) 
(dolist (note2 (elt ocm-chord-list (1+ i))) 
(if (and (or (= (third note2) 
(third notel)) 
(= (third note2) 
(third note3))) 

C< o * 

(mod (- (second notel) (second note2)). 
12) 

(mod (- (second notel) (second note3)) 

. X > V 

VZJ J ) 

(unless (remove-if 
tt'null 
(mapcar 

#' (lambda (note) 

(/= (second note) 

(second note2))) 
(remove-if 
t'null 
(mapcar 

#' (lambda (note) 

(if (= (third note) 

(mod (- (third notel) 1) 
7)) 

u— *\ *\ 

(elt ocm-chord-list (1+ i)))))) 
(setf (third note2) 

(mod (- (third notel) 1) 7)))))))))) 

; upward passing notes 
(ocm-chord-list 
(do* (Ci 0 (1+ i))) 

((= i (- number-of -chords 2)) 
ocm-chord-list) 
(dolist (notel (elt ocm-chord-list i)) 

(dolist (note3 (elt ocm-chord-list (+ i 2))) 

(if (= (third note3) (mod (+ (third note!) 2) 7)) 
(dolist (note2 (elt ocm-chord-list (1+ i))) 
(if (and (or (= (third note2) 
(third notel)) 
(= (third note2) 
(third note3))) 

(< 0 

(mod (- (second note3) (second note2)) 
12) 
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(mod (- (second note3) (second notel)) 
12))) 
(unless (remove- if 
#'null 
(mapcar 

#' (lambda (note) 

(/= (second note) 

(second note2))) 
(remove-if 
#'mill 

(mapcar # ? (lambda (note) 

(if (= (third note) 

(mod (+ (third notel) 
1) 

7)') 

note)) 

(elt ocm-chord-list (1+ i)))))) 
(setf (third note2) 

(mod (+ (third notel) 1) 7)))))))))) 
(morph-list (mapcar #'third (apply #'append ocm-chord-list))) 
(morphetic-pitch-list 
(mapcar #' (lambda (chromatic-pitch morph) 

(let* ((morphetic-octavel (floor chromatic-pitch 12)) 
(morphetic-octave2 (+ 1 morphetic-octavel)) 
(morphetic-octave3 (- morphetic-octavel 1)) 
(mpl (+ morphetic-octavel (/ morph 7))) 
(mp2 (+ morphe tic-octave 2 (/ morph 7))) 
(mp3 (+ morphetic-octave3 (/ morph 7))) 
(chroma (mod chromatic-pitch 12)) 
(cp (+ morphetic-octavel (/ chroma 12))) 
(difference-list (list (abs (- cp mpl)) 
■c (abs (- cp mp2)) 

(abs (- cp mp3)))) 
(morphetic-octave-list (list morphetic-octavel 

morphetic- octave 2 
morphetic-octave3) ) 

(best-morphetic-octave 
(elt morphetic-octave-list 

(position (apply #'min difference-list) 
difference-list)))) 
(+ (* 7 best-morphetic-octave) morph))) 
chromatic -pitch- list 
morph-list)) 

(opd (mapcar #' (lambda (tpcd-triple morphetic-pitch) 
(list (first tpcd-triple) 

(list (second tpcd-triple) 

morphetic-pitch) 
(third tpcd-triple))) 
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s ort ed- input -r epre sent at ion 
morphetic -pitch- list) ) 
(opnd (mapcar #' (lambda (opd-datapoint) 

" (list (first opd-datapoint) 

(p-pn (second opd-datapoint)) 
(third opd-datapoint))) 

opd))) 

opnd) ) 

(defun chromatic-pitch-chroma (chromatic-pitch) 
(mod chromatic-pitch 12)) 

(defun vector-less-than (vl v2) 
(cond ((null v2) nil) 
((null vi) t) 

((< (first vl) (first v2) ) t) 

((> (first vl) (first v2)) nil) 

(t (vector-less-than (cdr vl) (car v2))))) 

(defun p-pn (p) 

(let* ((m (p-m p)) 

(1 (elt '("A" "B" "C" "D n "E" "F" "G") m)) 

(gc (p-gc p)) 

(cdash (elt '(0 23 57 8 10) i)) 
(e (- gc cdash)) 
(i ) 

(i (cond ((< e 0) 

(dotimes (j (- e) i) 

(setf i (concatenate 'string i "f")))) 
((> e 0) 
(dotimes (j e i) 

(setf i (concatenate 3 string i "s")))) 

((= e 0) "n"))) 
(om (p-om p)) 

(oasa (if (or (= m 0) (= m D) 
om 

(+ 1 om))) 
(o (format nil ""D" oasa))) 
(concatenate 9 string 1 i o) ) ) 

(defun p-om (p) 
(div (p-pm p) 7)) 

(defun p-pm (p) 
(second p)) 



(defun div (x y) 
(int (/ x y))) 



(defun int (x) 

(values (floor x))) 

(defun p-gc (p) 
(- (p-pc p) 

(* 12 (p-om p)))) 

(defun p-pc (p) 
(first p)) 

(defun p-m (p) 

(bmod (p-pm p) 7)) 

(defun bniod (x y) 
(- x 

(* y 

(int (/ x y))))) 
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CLAIMS 

A method for computing the pitch names (i.e., Cf4, Bb3, etc.) of notes in a rep- 
resentation of music in which at least the onset time and MIDI note number (or 
chromatic pitch) of each note is given, comprising the steps of 

(a) computing for each pitch class 0 < c < 11 and each note n in the input, the 
pitch letter name S(c, n) G {A,B,C, D, E, F , G}, that n would have if c were 
the tonic at the point in the piece where n occurs (assuming that the notes are 
spelt as they are in the harmonic chromatic scale on c); 

(b) computing for each note n in the input and each pitch class 0 < c < 11 
a value CNT(c,n) giving the number of times that c occurs within a context 
.surrounding n that includes n, some specified number K p „ of notes immediately 
preceding n and some specified number if post of notes immediately following 



n; 



(c) computing for each note n and each letter name I, the set of chromas C(n, I) - 
{c | 5(c, n) = 1} (that is, the set of tonic chromas that would lead to n having 
the letter name t); 

(d) computing N(l,n) =. £ CNT(c,n) for each note n and each pitch letter 

ceC(n,l) 

name I; 

(e) computing for each note n, the letter name / for which N(l,n) is a maximum. 

The method of Claim 1, adapted to correct those errors in the output of the method 
of Claim 1 in which a neighbour note or passing note is erroneously predicted to 
have the same letter name as either the note preceding it or the note foUowing it, 
comprising the further steps of 

(a) lowering the letter name of every lower neighbour note for which the letter 
name predicted by the method of Claim 1 is the same as. that of the preceding 
note; 
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(b) raising the letter name of every upper neighbour note for which the letter name 
predicted by the method of Claim 1 is the same as that of the preceding note; 

(c) lowering the letter name of every descending passing note for which the letter 
name predicted by the method of Claim 1 is the same as that of the preceding 

. note; 

(d) raising the letter name of every descending passing note for which the letter 
name predicted by the method of Claim 1 is the same as that of the following 
note; 

(e) lowering the letter name of every ascending passing note for which the letter 
name predicted by the method of Claim 1 is the same as that of the following 
note; 

(f ) raising the letter name of every ascending passing note for which the letter 
name predicted by the method of Claim 1 is the same as that of the preceding 
note. 

3. Computer software or hardware adapted for performing the method of any preceding 
Claim 1-2. 
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ABSTRACT 

METHOD OF COMPUTING THE PITCH NAMES 
OF NOTES IN MIDI-LIKE MUSIC 
REPRESENTATIONS 

The invention described here consists of an algorithmic method called pslS that reli- 
ably computes the correct pitch names (e.g., C«4, Bb5 etc.) of the notes in a passage of 

„. j.! j n /ttt^t v,„mW of oaf.1i nnt.e in the 

tonal music, wnen given oniy me onset-bime w^j- — 

passage. 

The psl3 algorithm has been shown to be more reliable than previous algorithms, 
correctly predicting the pitch names of 99.81% of the notes in a test corpus containing 
41544 notes and consisting of all the pieces in the first book of J. S. Bach's Das Wohltem- 
perirte Klavier (i.e., pslS incorrectly predicted the pitch names of only 81 notes in this 
test corpus). 

Three previous algorithms (those of Cambouropoulos (1996, 1998, 2000, 2001, 2002), 
Longuet-Higgins (1976, 1987, 1993) and Temperley (1997, 2001)) were run on the same 
corpus of 41544 notes. On this corpus, Cambouropoulos's algorithm made 2599 mistakes, 
Longuet-Higgins's algorithm made 265 mistakes and Temperley's algorithm made 122 
mistakes. As pslS made only 81 mistakes on the same corpus, this provides evidence in 
support of the claim that pslS is more reliable than previous algorithms that attempt to 
perform the same task. 
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Figure 1: Examples of notes with identical MIDI note numbers being spelt differently in 
different tonal contexts (from Piston, 1978, p. 8). 



THIS PAGE BLANK (USPTO) 



THIS PAGE BLANK'ojsptoj 



< 7. 2 : D4- 25 U,Wm* : (6,-F4, 2), . (8; G4„ 2);. (10, E4, 2); f 2, F4, 1) 

(30, ■ W ■ -pf^- gTl'S" S'H'S' 

44 -B4 -15 (45, A4EE):,- : (46:V:GS4, 1), (:47,.rA4,.. 1), ; (26, A3, 2),. ^(28,-^:2 , 

fo^iM', (slv^h 'mmU »C4 V .VV, ^37,A3,1), :(38v G S3,^ 
(39; A3,1), (40, E4;,4), (4-4-B4-, 4) ■ '}- 

Eig U re 3: (a):' Ba TS : 1 to 4 of Bath's': Fugue- in D minor tor Root 1 of Das- WMllemperirfe 
Klavier (B:WV,g51i).. (b) The :<9PM? representation of tbe-scoxe in (a). 



•v. 
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1 S<-() 

2 Q «- (0, 1, 1, 2, 2, 3, 3, 4, 5, 5, 6, 6) 

3 for ct- 0 to 11 

4 mo <— 0 

5 s«- () 
6 

7 f or i <- 0 to |C| - 1 

8 s^-s ®iQ[(0\i] ^c)mod:l2]>,,:\ 

9 m 0 <- s[0] 

10 for i «- 0 to |C| - 1 

11 s' <- 0 <(s[i] - mo) mod 7) 

12 5^-50 («') 

13 return 5 



Figure 4: Algorithm for computing the spelling table 5. 
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1 "*«-<) 

2 ->• f or i <— 0 to.. |C| — 1 

3 ■ V<- (0,0, 0,0,0,0,0) 

4 

5 f or fc «- 0 to 11 

6 ti^fJL®(S[k}[%[) 

7 for j *— 0 to 11 

8 Sprev < 

s «-*p«r + *Mb'] 

10 R<— R® (PDS(max(y) , V)) 

11 return R 



Figure 5: Algorithm for computing the relative morph 
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1 H <- «O[0]» 

2 for i *- 1 to \C\ - 1 

3 if O[i][0]=O[i-l][0] ;;..v ( :\ 

4 ff[|if| 3 .l] *- H[|fl-| - i] ^W)> 

5 ©Is© 

6 " H <r- H @ ((0[i])) 

7 return i? 



Figure 6: Algorithm for computing the chord list 
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Figure 7: Neighbour note and passing note errors corrected by pslS. 
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1 f or i <- 0 to \H\ - 3 

2 ■ C-0 

3 for k <-0 to |fft» + 2]| - 1 

4 C-C©<-Hl* + 2]t*][1.3)> 

5 for ni ^ 0, to- •• •• 

6 if JT(s|[ni][l,3] € C 

7 for n 2 «- 0 to |fl"[t + 1]| - 1 

8 if H[i + 1] [na] [2] = [n x ] [2] 

9 if (i?[t + l][n 2 ][l]--H'[i][ni][l])modl2€{l,2} 

10 + 1) [n 2 ] [2] «- (ff [i + 1] [n 2 ] [2] + 1) mod 7 

11 if {H [i][nx][l] - H[i + l][n 2 ][l]) mod 12 e {1, 2} 

12 J=r[i + 1] [n 2 ] [2] <- (£T[t + 1] [n 2 ] [2] - 1) mod 7 

13 return i? 



Figure 8: Algorithm for correcting neighbour note errors. 
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V for i <- 0 to \H\ - 3 

2 for nj f- 0 to |i3"[i]| - 1 

3 for n 3 <-0 to |ii"[i + 2]| - 1 

4 if H[i + 2)[n 3 ][2] = (ff[*][m][2] - 2) mod 7 

5 f or U2 «- 0 to \H[i + 1]| - 1 r „ „. 
6- if H\i + l][n 2 }[2]= jr[i]K][2} or H[i + l][n 2 }[2] = H[i + 2}[n 3 ][2) 

7 if 0 < (flfflUlfl] - H[i + UtnalW) n»od 12 < - JT[i + 2] MM) mod 12 

8 £-<> 

9 for j *- 0 to \H[i + 1]| - 1 

10 if H[i + l]\j][2] = (fllilWPl - 1) mod 7 

11 £^<;®(H[i + l]\j]) 

12 0^-0 ,-v. 

13 for j *-0 to 1 '..,>■ - ■'- " 

14 if i?[i + l][n 2 ][l]/C[i]W 

is 0^0®«L?1> 

16 if 0 = 0 

17 + l][n 2 ][2] - (flMnJpJ - 1) mod 7 

18 return if 



Figure 9: Algorithm for correcting descending passing note errors. 
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( 1 f or i «- 0 to \H\ - 3 . 

' 2 for n x i- 0 to \H[i}\ - 1 

3 for n 3 <- 0 to \H[i + 2]| - 1 

4 if fi"fi + 2] [n 3 ] [2] = (H\x\ [m] [2] + 2) mod 7 

5 for n 2 <-0 to + 

6 if + 1] [n 2 ][2] = [n!][2] or H{i + 1] [n 2 ]{2] = H[i + 2] [n 3 ] [2] 

7 if 0 < (Jff[i + 2][n 3 ][l] - H[i + l][n 2 ][l]) mod 12 < (iTfi + 2][n 3 ][l] - H[z][m][l]) mod 12 

8 C«-<> . .. ... ■: :■ - - • 

9 for j <— o to |/3-[i + i]| - 1 .; ,. : 

10 if J ff[i + l][j][2] = (ff[i]M[2]+l)mod7 

11 C*-C©<-H"[i + lPl> 

12 0^-<> 

13 for j <-0 to ICI-1 

14 if J?[i+iJM(i]#Cb][i] 
is 0 0 e (CLfl) 

16 if 9 = 0 

17 -fffi + 1] M [2] «- (if [i] [n x ] [2] + 1) mod 7 

18 return if 



Figure 10: Algorithm for correcting ascending passing note errors. 
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1 <7 WO 

2 for i <— 0 to - 1 

3 . «- o' e ri[i) .... 

4 M'^-() 

5 for i <- 0 to |C| - 1 

6 M' *- M'e (0'[*][2]> 

7 return M' 



Figure 11: Algorithm for computing M' 
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1 p«-<) 

2 for i <- 0 to | J| - 1 

3 0l LJ[i][l]/12j 

4 o 2 <— 1 + Oi 

5 o 3 <— Oi — 1 

6 pi <- O! + M'[i]/7 
.7 p 2 «- o 2 + M'[i]/7 

8 . p 3 t-o 3 + M'[i]/7 

9 ' c <- J[i*i(l] mod 12,. 

10 j/ <- Ox + c/12 

11 I? -f— <|p' -Pi|, |p' - Pa!-, |p' - Ps!) 

12 w <— (o!,o 2 , o 3 ) 

..o - , vrDnoUinCnl DM 

14 Pf-Pe <M'[i] + 7o> 

15 return P 



Figure 12: Algorithm for computing P. 



THIS PAGE BLANK wro) 



13/13 



PPN(p) 


1 


m 4— mod 7 


2 


L <- ("A", "S" , "C", U U 


3 


l^L[m) 


4 


jc^-pM-^CbWAj) 


5 


A <- (0,2,3,5,7,8, 10). 


6 


c 7 A[m] 


7 


e *- £ c — C 


8 


z < — 


9 


if e < 0 


10 


for j <— 0 to • — e — 1 


11 




12 


else 


13 


if e> 0 


14 


for j<-0 to e- 


15 


i <- i © "s" 


16 


else 


■4 T 

-L / 




18 


o m - b[l]/7j 


19 


if m = 0 or m = 1 


20 


o <— o m 


21 


else 


22 


o <— 1 + o m 


23 


o str STH(o) 


24 


return Z © i © Ostr 



Figure 13: PPN algorithm. 
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